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We present a simple method to determine the mutual inductance M between two coils in a coupled 
AC circuit by using a digital dual-phase lock-in amplifier. The frequency dependence of the real 
and imaginary parts is measured as the coupling constant is changed. The mutual inductance M 
decreases as the distance d between the centers of coils is increased. We show that the coupling 
constant is proportional to d~" with an exponent n (« 3). This coupling is similar to that of two 
magnetic moments coupled through a dipole-dipole interaction. 
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I. INTRODUCTION 



II. BACKGROUND 



Faraday's law of magnetic induction states that a 
changing magnetic flux through a coil of wire with re- 
spect to time will induce an EMF in the wire. In a 
coupled circuit, there are two coils of wire; a primary 
coil that is connected in series with the voltage source, 
and a secondary coil that is not connected to any voltage 
source. The secondary coil receives energy only by induc- 
tion. The EMF in the secondary coil affects the voltage 
across the primary coil due to reflected impedance.'i'^l 



Here we present a simple method to measure the fre- 
quency dependence of the real and imaginary parts of the 
output voltage across a series resistance in the primary 
circuit of the coupled coils. Changing the separation dis- 
tance between the two coils facing each other leads to the 
change in the mutual inductance. The present method 
allows one to obtain a great deal of data in a reasonably 
short period of time. A background for the AC analy- 
sis of the AC coupled circuit in the frequency domain 
is presented in Sec. |lT] The frequency dependence of the 
real and imaginary parts of output voltage is formulated 
and is simulated using Mathematica. Our experimental 
results are reported in Sec. |IV[ The mutual inductance 
is determined as a function of the distance between the 
centers of two coils from the frequency dependence of the 
real and imaginary parts. We show that the mutual in- 
ductance changes with the distance d as d"" where n « 3. 
This coupling is the same as that of magnetic moments 
which are coupled by a dipole-dipole interaction. 



The detection technique of the signal in the present 
method is similar to that used in the nuclear magnetic 
resonance,'^ ^ where the behavior of the real part (the 
dispersion) and imaginary part (absorption) of the Bloch 
magnetic susceptibility can be measured as a function of 
frequency near the resonance frequency. 



A. Self inductance and mutual inductance of coils 
used in the present work 

Figure [ija) shows a photograph of coils which are used 
in the present work. For simplicity, we assume that each 
coil has a cylindrical form with a radius Rav and a length 
I. The separation between the centers of two coils (see 
Fig. [l]^a)) is the distance d. The number of turns of 
the coils 1 and 2 is the same {N). The magnetic fleld 
produced by coil 1 (current Ii) at the center of coil 2 is 
given by 
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for d ^ Rav using the Bio-Savart law.l^^l The voltage 
induced in the coil 2 is 
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where the magnetic flux is $ = BA and A is the cross- 
sectional area of the solenoid: A — ttR^^. Thus the 
mutual inductance M defined by V2 = —Mdli/dt is given 

by 
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The ideal self-inductance Lq is given by 
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for sufficiently long I, where I is the length of solenoid, 
/Xo (= 47r X 10"'' Tm/A) is a permeability, and N is the 
total number of turn. Thus we have 
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since M = ULq (see Eq.( 14) for general definition of M). 
The constant k is dependent only on the geometry of the 
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FIG. 1; (Color online) (a) A picture of two coils used in 
the measurement. Each coil (Heath Company Part No. 40- 
694) has a cylindrical form with the inner diameter (9.30 
cm), the outer diameter (12.5 cm), and the length (9.0 cm). 
The number of turns of the coil is A'^ = 3400. The value of 
Li = L2 are determined from the resonance frequency with 
Ci — C2 ~ 0.0038 ^F. The resonance frequency /o is equal 
to fo = l/(27rVIC) = 2850 Hz. An AC voltage was sup- 
plied by the lock-in amplifier, (b) Frequency domain diagram 
of AC coupled circuit used in the present experiment. The 
primary circuit is on the left side, the secondary circuit on 
the right side. The coils 1 and 2 are used in place of the in- 
ductors Li and L2. Li = L2 = 0.821 H. Ci = C2 = 0.0038 
fiF. Rio — R20 = 31 fl. RiL and R2L are resistances of coils 
1 and 2. Ril = R2L = 66 ft. Ri = Rio + Ril = 97 Q. 
R2 = -R20 + R2L = 97 n. 



coils. Note that our coils used in the present measure- 
ment have N — 3400 turn and I =9.0 cm. When Rav 
— 5.45 cm, we have = 133.7, where d is in the 

units of cm (see FigQa)). The self-inductance Lq can be 
calculated as 



I 



= 1.506 H. 



This value of Lq is larger than the actual value of L^xp 
(= 0.821 H, which will determined experimentally). The 
difference between L^xp and Lq is due to the deviation 
of the system from ideal one because of the finite length 
of the coil: L^xp = K^Lq. Km is called the Nagaoka 
coefficient and is defined aJ^ 
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R and I are the radius and length of coil, respectively, 
and K{a) and £'(0;) are the complete elliptic integral of 
the first and second kind. 
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The value of Km can be calculated as a function of 2R/1. 
The ratio Km experimentally obtained in the present 
work is L^xp/Lq = 0.821 H/1.506 H = 0.5458. This ra- 
tio is close to the ratio Km {— 0.6456) calculated from 
Eq.(|6]) for the ratio 2Rav/l = 2 x 5.45/9.0 « 1.211. 



B. AC coupled circuit 

In Fig. [ijb) we show the AC coupled circuit in the fre- 
quency domain.El The currents and voltages are all com- 
plex numbers. Using Kirchhoff's law, we can write down 
two equations. 
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where j = \/— 1, lu (= 27r/) is the angular frequency, 
Ii and I2 are the loop currents of the primary and sec- 
ondary circuit, E is the source voltage, Zi and Z2 are 
the impedance of the primary and secondary circuits, re- 
spectively. 



Zi — Ri + jXi,Xi — bjLi - 

Z2 = R2+ jX2,X2 = UJL2 
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and M is dependent on the distance between Li and L2 
coils. From Eqs.(|8]) and ([9| we have 



E = Z'iIi,Z[ — Z\ 



(10) 
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where Z'^ is the effective impedance of the primary cir- 
cuit. The effective impedance is rewritten as 



For simplicity, we assume the symmetric configuration 
such that Ri = R2 — R, Ci = €2 = C, and Li = L2 = L. 
Then we have X — Xi = X2 = ujL—I/ujC. The effective 
impedance Z[ can be written as 



Z[ = R{1 



)■ 



The voltage across Rio between AG in Fig.[l|a), Vr, is 



Vr — hRio — E—^G, 
K 



where 
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We define the ratio x = lj/ljo, where lvq = 1/ y/LC and 
X is always positive. The quality factor of the circuit is 
given by 
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The mutual inductance M is related to the self induc- 
tance L by 



M = ky/LiL2 = kL, 



(14) 



where fc is a constant and is smaller than 1. By us- 
ing these relation, C and L can be expressed by C = 
l/{ujoQR) and L = QR/ luq. Then G can be rewritten as 



G = (1 + jv, 
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which depends only on x, Q, and k. According to 
Eq.(17l, V becomes zero when x = and 1. The other 



possible x's for giving v — Q ca,he examined by the fol- 
lowing quadratic equation: 



Q2(i _ e)x^ + 2g2)x2 + Q2 ^ 0. 
The solution of this equation is formally given by 
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The values of ^ and v depend on x and are described as 
coordinates in the (/i, v) plane for convenience; (/x, u) — 
(0,0) at X = and 00, (1/2,0) at x = xi, and (1/(1 + 
fc2Q2)^0) at X = 1. 

We now consider only the case of fc < 1 and Q ^ 1, 
which corresponds to the present experiment. 

(1) State-I: 4g2(fc2g2 _ 1) _^ ^ > q . This condition 
is nearly equivalent to kQ > 1 since g ^ 1. There are 
two solutions xi and X2 (0< X2 < 1 < xi) besides x = 
and 1. In Fig. [2j^a) we show the simulation plot of the 
trajectory of the point (/i,i^) for kQ > 1 and Li ~ L2 
(the symmetric configuration), when x varies from x = 
to 00. This figure corresponds to the case of kQ — 1.515 
where Q = 151. 5 and fc = 0.01. 

(2) State-II: 4g2(fc2Q2 _ 1) + 1 < . This condition is 
nearly equivalent to kQ <1. There is no solution, besides 
X = and 1. In Fig. [2f^b) we show the simulation plot of 
the trajectory of the point (/i, v) for kQ < 1 and Li — L2, 
when X varies from x = to 00. This figure corresponds 
to the case of kQ = 0.758, where Q = 151.5 and fc = 
0.005. 

Figures |3ja) and (b) show simulated plots of the real 
part (/z) and the imaginary part (i/) as a function of x, 
where Q = 151.5 and Li = L2. The coupling constant 
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FIG. 2: (Color online) (a) State-I. Simulation plot of the tra- 
jectory denoted by the point (/i, v) for kQ > 1 and Li = L2 
(the symmetric configuration), when x {— f / fo) varies from 
a; = to 00. The figure corresponds to the case of kQ = 1.515 
where Q = 151.5 and k = 0.01. The point on the trajectory is 
located at (fi, v) = (0, 0) for 2; = 0, at (1/2, 0) for x = X2, at 
(1/[1 + [kQf], 0) for a; = 1, at (1/2, 0) for x = 2:1, and at (0, 
0) for a; = 00. (b) State-II. Simulation plot of the trajectory 
denoted by the point (/i, for kQ < 1, when x varies from x 
= to 00. The figure corresponds to the case of kQ = 0.7575 
where Q = 151.5 and k = 0.005. The point is located at (/i, v) 
= (0, 0) for s = 0, at (1/[1 + [kQf], 0) for x ^ 1, and at (0, 
0) for x = 00. 



k is changed as a parameter: fc = - 0.08. In Fig. |3ja), 
the double peaks in the ^ vs x curve are symmetric with 
respect to x = 1 and become closer and closer together as 
k is decreased and become a single peak at kQ « 0.5757. 
The imaginary part v (see Fig.jsf^b)) has a positive local 
minimum at a; = 0.9967 and a negative local maximum 
at X = 1.0033 in the limit of k ^ 0. Figure|3]^c) shows the 
trajectory denoted by the point (fJ-jiy) for Q — 151.5 and 
Li = L2 (the symmetric configuration), when x is varied 
from X = to c». The coupling constant k is changed as 
a parameter: fc = - 0.1. There is a drastic change of the 
trajectory from the state-I and state-II at kQ = 1, when 
kQ is decreased. For A: = 0, the trajectory is a circle of 
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FIG. 3: (Color online) Simulation plot of (a) the real part (/i) 
and (b) the imaginary part {u) as a function of x (calcula- 
tions), where Q = 151.5 and L\ = L2 (symmetrical case), k 
is changed as a parameter: fc = - 0.08. (c) Typical trajectory 
denoted by the point (/i, v) for Q = 151.5 and Li = 1/2, when 
X is varied from a; = to 00. is changed as a parameter: k 
= - 0.08. 



radius 1/2 centered at /i = 1/2 and v = Q. 
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III. EXPERIMENTAL PROCEDURE 

The AC coupled circuit in the frequency domain was 
configured in Fig.[TJb). An AC voltage was placed across 
the primary circuit, Vs{t) — Re[Ee^'^*], where E is the 
complex voltage source. Two LCR circuits are placed 
together with inductors, resistors and capacitors. In 
Fig. [l|a) we show the overview of two coils used in the 
present measurement. The minimum distance between 
the centers of the coils is 10.2 cm. The two inductors are 
placed on a meter long track. 

We used a digital dual-phase lock-in amplifier (Stan- 
ford, SR850),'^ which is set at both a dual-phase mode 
and a frequency-scan mode. The real part fj, (the in-phase 
signal), and imaginary part (the out-of-phase signal) 
were measured simultaneously, when the frequency was 
continuously changed across the desired frequency range. 
The frequency range initially chosen was 2500 Hz to 3500 
Hz. At the 25.2 cm separation the frequency scan was 
only needed to be between 2700 Hz and 3000 Hz. The 
measurements are repeated for numerous distances rang- 
ing from d = 10.2 to 70.2 cm. The output voltage across 
the resistance (AG) in Fig. [ijb) is given by 



Voutit) 



Re[{fi + jv)E'^ exp[j(wt + 0o)] 
EQiJ,cos{u!t + ipo) + E'i^vcos{u)t + i 



>o + vr/2), 
(22) 



where E = Eoe^"^", is the phase, and Eq = EqRiq/R. 
The in-phase component of the lock-in amplifier is equal 
to Eq^ and the out-of phase component is equal to EqV, 
where Eq — 4.97 mV. Consequently, one can determine 
the values of the real part fi and the imaginary part v 
defined by Eqs.(16l and (17 1, independently. 



IV. RESULTS 

We have measured the frequency dependence of the 
real part fi and the imaginary part v when the distance 
d between the centers of two coils is changed as a param- 
eter: d = 10.2 - 70.2 cm. Our results are shown in Figs.|4] 
- 16] Figures |4] and |5] show the experimental plots of the 
real part (/it) and the imaginary part (v) as a function 
of X (= f/fo), where /o = 2850 Hz and Q = 151.5. In 
Figs. |4|a) and (b) double peaks of /z become closer and 
closer when d is increased and the double peaks become 
a single peak around d — 36.2 cm. The double peaks are 
not symmetric with respect to a; = 1. The peak at the 
lower- a; side is higher than that at the higher-x side. The 
real part has a local minimum at x which is a little larger 
than 1. In Fig.|5]^a) and (b) the imaginary part ly crosses 
the V = line at a: = X2, x ~ 1, and x — xi {> X2). 
The positions xi and X2 become closer and closer as d is 
increased and combine into the position a; w 1, but not 
at X = 1. Figure [6] shows the experimental trajectories 
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FIG. 4: (Color online) Experimental plot of (a) the real part 
(/i) and (b) the imaginary part (v) as a function of x {— f/fo) 
(in the frequency scan). The distance d between the centers 
of two coils is changed as a parameter, d = 10.2 cm - 17.2 
cm. fo = 2850 Hz. Q = 151.5. 



of the point {^,v), when x (= f/fo) is varied from x = 
0.947 (/ = 2700 Hz) to 1.053 (/ = 3000 Hz). The dis- 
tance d is changed: d — 13.2 — 36.2 cm. The transition 
occurs between the state-I and state-II at d « 27.2 cm. 
The overview of our trajectory is similar to the simula- 
tion plot as shown in Fig. ^c). However our trajectory 
rotates clockwise compared to the ideal simulation plot 
(the symmetric configuration). The deviation of our tra- 
jectories from the ideal case (Fig. |3jc)) is partly due to 
the asymmetric configuration (Li is slightly larger than 

L2)- 

Figure [7]^a) shows the zero-crossing frequencies nor- 
malized by fo for the imaginary part (v) (at which v be- 
comes zero) as a function of the distance d (cm), where 
X2 = /2//0 (< 1), Xc = fc/fo (« 1), and xi = /1//0 
(> 1). The value of little different from 1. In 

Fig. ^a) we also show the normalized frequency defined 
by (a;ia:2)^''^. This frequency decreases with increasing 
d. This implies that the parameter k decreases with in- 



creasing d as predicted from Eq.(21 1 



Figure ^h) shows the normalized peak frequencies of 
the real part (/x) as a function of the distance d (cm), 
where Xd = fd/fo (< 1) and a;„ = fu/fo (> !)■ The 
real part takes double peaks at the lower and upper 
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FIG. 5: (Color online) Experimental plot of (a) the real part 
{fi) and (b) the imaginary part (i^) as a function of x {— ///o) 
(in the frequency scan). The distance d between the centers 
of two coils is changed as a parameter, d = 19.2 cm - 36.2 
cm. fo = 2850 Hz. Q = 151.5. 
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FIG. 6: (Color online) Experimental trajectories of the point 
for Q = 151.5, when x (= f/fo) is varied from x = 
0.947 (/ = 2700 Hz) to 1.053 (/ = 3000 Hz), fo = 2850 Hz. 
The distance is changed as a parameter: d = 13.2 - 36.2 cm. 
The deviation of the experimental trajectories from the ideal 
one as shown in Fig. [sj^c) is partly due to the asymmetric 
configuration (Li is slightly larger than L2). 



FIG. 7: (Color online) The normalized zero-crossing frequen- 
cies of the imaginary part (u) (at which v becomes zero) as a 
function of the distance d (cm). X2 = ,f2//o (< 1), Xc = /c//o 
{k: 1), xi = /i//o (> 1), and fo = 2850 Hz. The normalized 
frequency defined by {xiX2y''^ is also shown for comparison, 
(b) The normalized peak frequencies of the real part (/i) as 
a function of the distance d (cm). Xd — fd/fo (< 1) and 
Xu = fu/fo (> 1). fo ~ 2850 Hz. The real part /i takes 
two peaks at the lower and upper frequencies fd and for 
d < 30 cm. The normalized frequency defined by (xuXd)^'''^ 
is also shown for comparison. 



frequencies fd and /„ for d < 30 cm. The d dependence of 
x„ and Xd is similar to that of xi and X2, respectively. In 
Fig. [7]jb) we also show the normalized frequency defined 
by {xuXdY^"^ as a function of d. This frequency decreases 
with increasing d like {xiX2)^^'^ in Fig. [Tia). 



V. DISCUSSION 

First we show that, from a view point of physics, the 
AC coupled circuit with the mutual inductance is equiv- 
alent to the magnetic moments Mi (= hqNIiA) for the 
coil 1 and A/2 {— fioNl2A) for the coil 2. They are cou- 
pled with a dipole-dipole interaction defined by 

1 , Mi-M2 3(Mi-r)(M2-r) , 

(^12-;; [ ^ E J, (^3) 

47r/io J' 

where r is the position vector connecting the centers of 
coils 1 and 2. When both Mi and M2 are parallel to 
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FIG. 8: (Color online) (a) and (b) Plot of fc as a function of 
the distance d. Q = 151.5. (a) The value of k is derived from 
the prediction that the real part (/i) is equal to l/[l + (fcQ)^] at 
X — 1 for the symmetrical configuration (Li = L2). The best 
fitted curve to the expression given by Eq.( |25[ ) is denoted by a 
solid line, (b) The values of k are derived from the prediction 
that xi and X2 are described by Eqs.( 19 1 and ( 20 1. The values 



of k are numerically solved for each d. The best fitted curves 
are shown by the dotted and solid lines in the figure. 



the direction of r, a parallel alignment of two magnetic 
moments is energetically favorable, 



Ul2=- 



2 M1M2 



(24) 



From the definition of the mutual inductance M, the in- 
teraction energy U12 can be described by U12 = — M/1/2, 
leading to the mutual inductance which is the same as 
Eq.(|3| derived from Faraday's law. 

It is predicted that the parameter k changes with dis- 
tance according to Eq.([5|; k is proportional to . The 
first method to determine the parameter fc as a function 
of the distance d, is as follows. As shown in Sec. |ll] 
it is predicted that in the symmetrical configuration 
(Li = L2 = L), the real part takes a ^ = + (f^Q)^] 
ai X — I. Note that the imaginary part v is equal to zero 
at a; = 1. Experimentally we determine the value of /i at 
x = 1 as a function of d. The values of k are derived from 
the above expression with Q = 151.5. Figure [sFa) show 



the plot of fc vs d thus obtained. The value of k drasti- 
cally decreases with increasing d and almost reduces to 
zero at d = 35 cm. The least-squares fit of the data of fc 
vs d to an expression 



fc = 



C 



(25) 



yields a constant ( = 38.5 ± 5.0 and the exponent n = 
2.60 ± 0.05, where d is in the units of cm. The value 
of C is rather different from the predicted value for the 
present coils = 151.5), while the value of n is rather 
close to the predicted value {n ~ 3). The large deviation 
of the experimental value of C from our prediction may 
be related to the asymmetric configuration of Li and L2 
in the present system, where L2 is slightly lower than Li 
(which will be discussed later). As shown in Fig. |4]^a), 
the value of x where the real part fi has a local minimum 
is not equal to a; = 1, and shifts to the high- a; side. 

The second method to determine the value of fc is as 
follows. In Sec.|IT] it is predicted that the imaginary part 
v takes zero-crossing at x — X2, 1, and xi in the case of 
the symmetrical configuration (ii — L2 — L). Note 
that the imaginary part v is not always equal to zero 
at X = I partly because of the asymmetric configuration 
in the present experiment. The value of fc for each d 
is derived by applying the Mathematica program called 
"FindRoot" to Eq.(19l with the experimental value of 
X2 and to Eq.(20l with the experimental value of xi (see 



Fig.|7][a)), since Eqs.(19l and (20) are complicated func 
tions of fc. In Fig. jsj^a) we show the value of fc as a 
function of d thus obtained. The value of fc drastically de- 
creases with increasing d. The value of fc is a little larger 
than those obtained from the first method at the same 
d. The least-squares fit of the data of fc vs d to Eq. ( 25 1 
yields the parameters C = 57.1 ± 6.9 and n = 2.48 ± 0.05 
for xi and C = 149 ± 33 and n = 2.92 ± 0.09 for X2. The 
latter result is in excellent agreement with the prediction 
(C — 151.5 and n — 3.0). Such different values of C are 
partly due to the effect of the asymmetric configuration of 
coils. Nevertheless, it may be concluded experimentally 
that two magnetic moments made from coils are coupled 
through the dipole-dipole interaction with the exponent 
n being equal to 3. 

Finally we discuss the effect of the asymmetric config- 
uration on the trajectory in the (fJ-jv) plane. As shown 
in Fig. [6] the trajectory rotates clockwise compared to 
the case of the trajectory in the symmetric configuration. 
Figure |9ja) shows the simulation plot of the fi as a func- 
tion of X for the asymmetric configuration (Li = 0.8201 
H and L2 — 0.8215 H) as fc is changed as a parameter. 
Double peaks of fi around a; = 1 are not symmetric with 
respect to a; = 1. The peak at the high- a; side is higher 
than that at the low- a; side. Double peaks become closer 
and closer as fc is decreased. Figure [9|b) shows the sim- 
ulation plot of the trajectory in the (/i,!') plane under 
the same condition as Fig. [9|a). The trajectory rotates 
counterclockwise compared to the case of the trajectory 
in the symmetric configuration. Figures 10 ^a) and (b) 
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FIG. 9: (Color online) (a) Simulation plot of the real part 
(/x) as a function of x, where Q = 151.5 and L\ — 0.8207 H 
and 1/2 = 0.8215 H (the asymmetric configuration), (b) The 
trajectory of the point {fJ,,v) for Q = 151.5, Li = 0.8207 H 
and L2 — 0.8215 H, when x is varied from a; = to 00. The 
coupling constant k is changes as a parameter: fc = - 0.08. 



FIG. 10: (Color online) (a) Simulation plot of the real part 
(/i) as a function of x, where Q = 151.5 and Li = 0.8207 H 
and 1/2 = 0.8198 H. The coupling constant k is changed as 
a parameter: k — - 0.08. (b) Typical trajectory denoted by 
the point {jJ.,1') for Q = 151.5, when x is varied from a; = 
to 00. The coupling constant k is changes as a parameter: k 
= - 0.08. 



VI. CONCLUSION 



show the simulation plot of /i as a function of x and the 
trajectory in the plane for the asymmetric config- 

uration (Li = 0.8201 H and L2 = 0.8198 H) where k is 
changed as a parameter. Double peaks of fi around x = 
1 are not symmetric with respect to a; = 1. The peak 
at the high- a; side is lower than that at the low- a; side. 
The trajectory rotates clockwise compared to the case of 
the trajectory in the symmetric configuration. These fea- 
tures are in good agreement with those observed in the 
present measurement (see Fig. [sjja) for the /i vs x curve 
and Fig. [6] for the trajectory). So we can conclude that 
Li is a little larger than L2, which means the asymmetric 
configuration for the present measurement. 



We present a simple method for determining the mu- 
tual inductance of the AC coupled circuit using a digital 
dual-phase lock-in amplifier. This method allows one to 
get a large amount of data on the frequency dependence 
of the real and imaginary part of the AC output volt- 
age in a reasonably short time. Our experimental results 
show that the coupling constant of the two coils is pro- 
portional to with an exponent n (« 3), where d is 
the distance between the centers of coils. This coupling is 
similar to that of two magnetic moments coupled through 
a dipole-dipole interaction. 
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